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We explore the formulation of non-rational 2D quantum gravity in terms of a chiral CFT on 
a Riemann surface associated with the target space. The CFT in question is constructed as 
the collective theory for a matrix chain, which is dual to a statistical height model on dynamical 
triangulations. The heights are associated with the sheets of the Riemann surface, which represents 
an infinite branched cover of the spectral plane. We consider two examples of height models: the 
SOS model and the semi-restricted SOS (SRSOS) model, in which the heights are restricted from 
below. Both models are described in the continuum limit by theories of 2D quantum gravity with 
conformal matter, perturbed by a thermal operator (1,3). We give a compact operator expression 
for the n-loop amplitudes as a collection of target space Feynman rules. The n-point functions of 
local fields are obtained by shrinking the loops. In particular, we show that the 4-point function of 
order operators in the SRSOS model coincides with the 4-point function of the "diagonal" world 
sheet CFT studied in Iffl. 



1. Introduction 

In a previous paper we studied non-rational theories of 2D Quantum Gravity (QG), 
having continuous spectrum of the matter central charge — oo < c < 1. We explored 
and generalized the ground ring approach to the construction of the tachyon correlation 
functions. The method was also tested for an unconventional model of 2D quantum gravity 
with interaction which induces matter charges corresponding to the diagonal of the Kac 
table. In this paper we develop a different technique for the computation of the tachyon 
correlators based on discrete realizations of the 2D gravity with non-rational conformal 
matter. 

The simplest statistical model that leads to a non-rational CFT is the SOS height 
model, which we will consider on a triangular lattice. The local fluctuating variable in 
the SOS model is an integer height x, which can jump by ±1 between neighboring sites. 
The Boltzmann weights are complex and depend on two parameters, the "background 
momentum" po G [0, 1] and the temperature coupling T > 0. By restricting the heights 
of the SOS model to positive integers we obtain another height model which we call semi- 
restricted SOS, or SRSOS model. The Boltzmann weights of SRSOS model are real and 
anti-symmetric with respect to the reflection x —* —x. As we will consider the two models 
alongside, we denote by X the target space of any of them, 

f Z for SOS; , . 

\Z+ for SRSOS. ^ ' 

At the rational points the SRSOS model can be further restricted to the RSOS models, 
which have positive Boltzmann weights The RSOS models form the A-series of the 

ADE Pasquier models In this sense the SRSOS model is "universal cover" of the 
A-series of minimal models. Its target space can be identified with the Dynkin diagram of 
the Aoo Lie algebra, considered as a certain limit of A^. It is also possible to construct the 
non-rational extension of the Dri-series, with symmetric weights, which we discuss briefly 
in Appendix C. A common property of the height models is that they can be mapped to a 
gas of self-avoiding and mutually avoiding loops on the dual lattice. The observables are 
introduced by assigning special weights to the non-contractible loops, while the activity of 
the contractible loops 

n 2 cos(7rpo) (1-2) 

is determined solely by the background momentum. The temperature coupling determines 
the "mass" of the loops. 

At the critical temperature T = Tc the large-distance behavior of the SOS model 
is that of a gaussian field with curvature term, while the SRSOS model is described by 
Virasoro CFT with charge reflection symmetry. Both CFT's have the same central charge 

c = 1 — 6 ^° (1.3) 

critical 1-|-J>q ' \ ' / 
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There is a set of primary fields whicfi can be defined at microscopic level. These are the 
order operators, which are realized by the eigenfunctions of the adjacency matrix of the 
target space X. Each order operator is characterized by a momentum p. The degenerate 
order operators, placed along the diagonal of the Kac table, have momenta p = mpo, 
m G Z_|_. The order operators are well defined also away from the critical temperature, 
where the long distance behavior of the height models is described by a perturbation with 
the thermal operator $1,3 

The height models coupled to gravity are defined by generalizing the Boltzmann 
weights, originally defined on a flat lattice, to an arbitrary triangulated surface, according 
to the prescription in P,p!0[|. Then the path integral over two-dimentional geometries is 
discretized as a sum over all triangulations with given topology. 

The SOS and SRSOS models have their dual formulations in terms of doubly infinite 



or semi infinite matrix chains, similar to the ADE matrix models fill. The collective field 



theory for the matrix chain is that of a chiral boson on a Riemann surface representing 
an infinite branch cover of the spectral plane of the random matrices. The sheets of the 
Riemann surface are labeled by the points x of the target space. The point on the Riemann 
surface are enumerated by pairs (z, x) where z G C is the projection on the spectral 
plane and a; G X labels the sheet. Thus each point of the Riemann surface corresponds 
to a boundary condition (FZZ-brane), with z being the complex boundary cosmological 
constant and x being the height at the boundary. At the rational points the Riemann 
surface represents the algebraic curve of the corresponding minimal model of 2D gravity 

We will be interested only in the genus zero amplitudes, which are well defined quan- 
tities even in a non-unitary theory. The disk and the annulus amplitudes are given respec- 
tively by the one-point and by the two-point functions of a free boson on the Riemann 
surface. The amplitudes with more than two boundaries are produced by the interaction 
terms. The latter, associated with the branch points, are needed to preserve the conformal 
invariance [|13[ . The operator solution of the Virasoro constraints yields a set of Feynman 



rules for calculating the loop correlation functions, which allows to evaluate any such func- 
tion by a finite sum over all possible intermediate channels. In the case of rational ADE 
string theories these rules were formulated in |n . 

We work out explicitly the example of the 4-point function and will compare the re- 
sult to the prediction of the world-sheet GET. In the case of the SOS model, our result 



reproduces at the critical point the Di Erancesco-Kutasov formula [|T5[ obtained for gaus- 
sian matter field with curvature term. In the case of the SRSOS model, our result for 
the 4-point function matches the one obtained in for the "diagonal" perturbation of 
Liouville gravity. This perturbation is generated by the Liouville-dressed vertex operator 
of dimension zero obtained from the identity by charge reflection. 
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The paper is organized as follows. In sect 2 we give the definition of SOS and SRSOS 
models on a triangulated surface with curvature defects and their representation in terms of 
a gas of loops. The loop gas representation of the correlation functions of order operators 
is sketched in Appendix A. In sect. 3 we reformulate the height models on dynamical 
triangulations in terms of doubly or semi infinite matrix chains. In sect. 4 we give the 
solution for the disk amplitude in the continuum limit. In sect. 5 we construct the target 
space CFT as a bosonic field theory on the Riemann surface associated with the disk 
amplitude. The Feynman rules for calculating the n-loop amplitudes are obtained in sect. 
6. In sect. 7 we work out the example of the 4-point function of order operators. Summary 
of the results is given in sect. 8. 

2. The SOS and SRSOS models on a triangulation with curvature defects 

2.1. Local Boltzmann weights and mapping to a loop gas 

Let ^ be a triangulated surface with the topology of a sphere with n boundaries. The 
triangulation is characterized by its vertices r, its links I =< rir2 > and its triangles 
A =< rir2r3 > which are assumed equilateral. The curvature defects are associated with 
the points r E Q with coordination number Cr (the number of triangles having this point 
as a vertex) different than 6. The local curvature at such point is equal to the deficit angle, 

^^ = ^(6-c^). (2.1) 
Similarly, the boundary curvature at the points r G dQ is 

i^r = |(3-c,). (2.2) 

By Euler's relation the total curvature is 

^Rr + 2^kr = 47r(2 - n). (2.3) 
reG reog 

The SOS and SRSOS models on the triangulation Q are defined as follows. To each 
node r G ^ we associate an integer height Xr G X, where the target space X is defined 
in ( p..l| ). The allowed height configurations are such that the heights of two nearest- 
neighbor points are either equal or differ by ±1. If one thinks of the target space as 
a one-dimensional graph (Fig. 1), then the allowed height configurations are the maps 
^ ^ X which preserve the nearest neighborliness, i.e., such that points are mapped to 
points and links are mapped either to points or to links. Since the space X is discrete, 
the momentum space P is compact: p ~ p -|- 2. In the case of SRSOS there is additional 
refiection symmetry, p ~ —p. 
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1 2 3 4 5 ... 

Fig. 1: The target space of the SRSOS model as a one-dimensional graph. 

The partition function Zg{xi, Xn) is the sum over all allowed maps ^ ^ X such that the 
boundaries have fixed heights xi, Xn- Besides the standard local factors Wa associated 
with the faces of Q, there are extra weights W, that come from the curvature defects: 

Zg{xi,...,Xn) = Yl WAiXr,, Xr^.Xr J Y[W,{Xr). (2.4) 

S— »X <rir2r3> r 

The Boltzmann weights are expressed in terms of the function 

= / exp(-i7rpoa;) for SOS; , . 

" \ sin(7rpox) for SRSOS ^ ' 

which plays the same role as Peron-Frobenius vector in the ADE height models. It is an 
eigenvector of the adjacency matrix of the graph X with eigenvalue 2cos(7rpo)- The back- 
ground momentum po G [0, 1] is assumed in this paper to be non-rational. By definition, 
the weight of a triangle A123 with heights xi,X2, xs is invariant under cyclic permutations 
and is non-zero only if the heights of each pair of points are either the same or adjacent. 
The possibilities are either xi = X2 = X3 or xi = X2 = X3 ± 1, up to a cyclic permutation, 
and have weights 

Wa{x,x,x) = 1; Wa{x,x,x±1) = ^(^^^ ' , (2.6) 

where T is a positive constant called temperature. The weight associated with the vertex 
r G ^ with curvature defects are 



Wix, 



(^.J^'-/^- if rG^; .2 7) 

(^,j^'^/2- ifr Ga^;' 



While Boltzmann weights of the SOS model are complex for any po 7^ 0, those of the 
SRSOS model are real, but not always positive. 

The height model is also described in terms of a loop gas. Let us represent the weights 
of the triangles with admissible heights symbolically as 



x-\ x+\ 



(2.8) 

That is, whenever the three heights round the triangle are not the same, a line separates 
the two vertices with the same height from the vertex with different height. The line is 
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made by two segments at angle 7r/3, orthogonal to the edges they start from. Then each 
admissible height configuration determines a pattern of closed nonintersecting loops (Fig. 
2). With our choice of boundary conditions the domain walls do not cross the boundary. 




Fig. 2: A loop configuration on a triangulation with boundaries. The boundary of 
the highlighted domain has rib = 3 components, two of them being dynamical loops. 

The Boltzmann weights in ( p.4|) are such that the weight of each loop gas configuration 
factorizes to a product of the weights of the connected domains of constant height. The 
weight of a domain T> is topological - it depends on the height x of the domain only 
through its Euler characteristics, i.e., the number of its boundaries^ 

Op(x) = = # boundaries of V. (2.9) 

The proof goes as follows. We can distribute the weights of the triangles that belong to two 
different domains, given by the second term of eqn. ( p.6|) , in such a way that the weight 
of each domain V is of the form {Sj^Y, where x is the height of the domain. It remains 
to evaluate the power q. It receives contributions from the vertices in the bulk inside the 
domain wall, from the vertices of the boundaries of Q (if any) that are also boundaries of 
and from the triangles along the domain walls that delimit T>. If there are no domain 
walls at all, that is, T> = then the total power of is g = 2 — n according to the Euler 
relation ( |2.3|) . If there are one or more domain walls, each domain wall contributes a factor 
S"^, where K is the total boundary curvature along it. Indeed, the weight of a domain wall 
is a product of factors {Sx)^^^^ that come from the triangles along it and account for the 
defects ±7r/6 along its edge. Applying again the Euler relation we find q = 2 — Ub, where 
Uf, is the total number of (connected) boundaries of the domain, given by the number of 
boundaries of Q trapped in the domain plus the number of the domain walls. 

^ In the case of a lattice with higher genus, there would be also a factor {Sx)^^ for each handle 
trapped in the domain. 
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As a consequence, the partition function ( p.4|) can be reformulated as a sum over 
all configurations of self and mutually avoiding loops on the dual lattice and heights x-r> 
associated with the connected domains V: 

= Yl ^"""^ E U^ryixjy). (2.10) 

loop configu- heights XT) 

rations on Q 

Here Ltot is the total length of the loops, which is equal to the total number of loop 
segments. The sum over the heights of the domains, {xx>}, is such that the heights of 
the neighboring domains differ by ±1. In this way the mapping to a loop gas, originally 
established for a fiat lattice Jl^, remains true also in the case of a lattice with curvature 



defects. 

For the disk partition function (n = 1) the one can perform readily the sum over the 
heights but one using repeatedly the relation 

2cos7rpo5'x- (2-11) 

x' 

In this way the sum over heights produces a factor (2cos7rpo)'^ loops ^Yme^ factor S^^ 
depending on the boundary height xi. Up to this last factor, this is the disk partition 
function of the 0(n) loop gas with activity n = 2 cosyrpo per loop. 

In the case of more than one boundary, the sum over the heights can be performed 
by Fourier transformation with the eigenvectors of the adjacency matrix 

(p)^re-P- for SOS; , . 

\ sin{Txpx) for SRSOS. ^ ' ' 

In this case the loops can have different activities. The contractible loops, i.e., the loops for 
which there is no topological obstacle for shrinking them to points, have activity 2 cosTrpo- 
A non-contractible loop contributes a factor (2cos7rp), where p G P is the momentum 
associated with the loop. In Appendix A we give the loop gas representation of the 
correlation function of n order operators on the sphere, which is the Fourier transform of 
the the partition function ( p.lO|) with boundaries of length zero. 

2.2. Phase diagram and critical points 

The critical thermodynamics of the SOS and SRSOS models is that of the 0(n) loop gas 
with n = 2 cos(7rpo)- The phase diagram (Fig. 3) of the 0(n) model on the the honeycomb 
lattice (dual to the regular triangulation) was first established in |]T^. At the critical 



temperature Tc = 2 cos the loop gas model is solvable and is described by a CFT with 
central charge 

^critical ^ "1 i ' (^2.13^ 

1 +Po 
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The critical point is also known as the dilute phase of the loop gas. For T > Tc the theory 
has a mass gap. The low-temperature, or dense, phase T < Tc is a flow to an attractive 
fixed point [jl^ at T^'^^^'^ = 2 sin , where the theory is again solvable and is described 
by a CFT with smaller central charge 



Ch =1-6 

cleiisc 



1 -Po 



(2.14) 



> > > — < < < < 
low temperature 
dense loops 



Oo 



-93- > 

high temp 
finite loops 



Fig. 3: Phase diagram of the loop gas on a regular (honeycomb) lattice. 

The scaling behavior of the model in the vicinity of a critical point is described by an 
action of the form [Bl 



A = Aritical + ST J $1,3 (2.15) 

where ST = T — Tc and $1,3 is the thermal operator with conformal dimensions 

1 -Po 



Ai,3 = Ai 



l+Po 



(2.16) 



When the thermal operator $1,3 is added to the action, it generates a tension of the loops 
proportional to 5T. For ST > the deformation (|2.15|) describes, going from short to the 
long distance scales, the flow to a massive theory with mass gap 



m 



(2.17) 



When 6T < the deformation ( |2.15| ) describes a massless flow between two different CFT 
with central charges ( p.l3| ) and ( |2.14| ). In the CFT for the dense phase, describing the IR 
limit, the flow to the attractive critical point is generated by another, irrelevant, operator, 

$3,1- 

2. 3. Mapping to a gaussian field 
• SOS model 

It has been argued that at the critical points the large distance behavior of the SOS 
model on a lattice with curvature defects is that of a gaussian fleld x(cri,a2) with a 



background charge proportional to the local curvature []16| , p!8[| . The product of the factors 
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( p77|) associated with the curvature defects yields the curvature-dependent term in the 
effective action 



^gauss 



^ J £a[{Vx? + ieox\- (2-18) 



The two critical points are described by different values of the background charge eo, 
related to the background momentum po by 

Co = ± =, with (+ for dilute, — for dense). (2 19) 

In the SOS model the order operators i(^p{x) = si^^ / Sr^ have the form of plane waves 

t{;p{x) = e*^(Po-P)^ = e^^*'^^. (2.20) 

At each of the critical points they can be identified with vertex operators of the gaussian 
field with conformal dimensions 

p'^ -pI q{q - Po) 



4(1 ±po) l±Po 



[+ for dilute, — for dense). (2-21) 



• SRSOS model 

In this case the correlation functions are constructed from vertex operators and screen- 
ing charges. To get an intuition about how the vertex operators appear in the microscopic 
theory, let us notice that in the SRSOS model the effect of the curvature does not reduce 
to a pure phase factor, as it was the case for the SOS model. Here we follow the argument 



of [|T8[. After formally expanding, with given by ( |2.5| ), one gets 

log^^ = -inpox - e^'^^Po^ + ... 

The exponential term corresponds to a vertex operator ( |2.2Q ) with q = po+ integer. The 
choice q = po + 1 gives one of the screening charges added to the action ( p.l8|) , e = a/1 -|- pq. 



3. SOS and SRSOS models on dynamical triangulations 

When considered on dynamical triangulations, the heights models give the desired solvable 
discrete models of non-rational 2D quantum gravity. The partition function in the ensem- 
ble of triangulations with the topology of a sphere with n boundaries, or shortly n-loop 
amplitude, is defined as 

^(Li,xi),(L2,x2),...,(L„,x„) = ^ Zg^^ ^^(Xi, . . . ,Xn). (3.1) 

The sum goes over all triangulations with the above topology and fixed lengths 

(number of edges) Li, . . . , L^^ of the boundaries, and the "cosmological constant" k coupled 
to the area A (the number of triangles) of the triangulation. 
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3.1. Dual description in terms of doubly infinite or semi infinite matrix chains 

There is a dual, matrix, formulation of the height models on dynamical triangulations in 
terms of a doubly infinite (for SOS) or semi-infinite (for SRSOS) matrix chains. These 
matrix models are similar to the ADE matrix models [|lll , with some subtleties related to 
the fact that they describe non-unitary theories. The fiuctuating variables are hermitian 
X -^x matrices M^; and the complex rectangular A^^; x A^^+i matrices (x G X). 
The hermitian matrices M^; are naturally associated with the sites x while the complex 
matrices C^; are associated with the links < x, x -|- 1 > of the target space. The partition 
function of the matrix chain is given by the integral 

Z[V] = j [dM] [dC\[dC^] exp ' (3-2) 

A=Y,Sx tr (iM^ -\Kml+ TClC - kC^CIM, - KM,+^ClC, \ (33) 

The perturbative expansion of the matrix integral is constructed according to the Feynman 
rules in Fig. 4, which correspond to the five terms in the action ( p.3| ). The correlation 
functions in the matrix model can be expressed in terms of fat (double lined) graphs. We 
will be interested in the limit 00 for all x G X, in which limit only the planar graphs 

survive. Each planar graph is in one-to-one correspondence with a triangulated surface. 
The vertices of the triangulation (dual to the faces of the planar graph) have label x G X. 
Thus the sum over planar diagrams reproduces the path integral for the height model in 
the ensemble of triangulations. 



x+l 



x+1 



x+l 



x+l 



Fig. 4: Feynman rules for the SRSOS matrix model 

In order to obtain the correct Boltzmann weights, we must tend the size of the matrices 
to infinity in a particular way: 



9s ^0. 



(3.4) 



Using the interpretation of the Feynman graphs as triangulated surfaces, one can easily see 
that the n-loop amplitude ( |3.1| ) is equal to the genus zero connected correlation function 
of a product of n traces 



(3.5) 
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The reader might object that the hmit (|3.4| ) does not exist, because is not a positive 
number. The hmit is in fact well defined only at the rational points of the SRSOS model, 
Po = l/{h + 1), with h G Z_|_. In this case Sx is Peron-Probenius vector for the adjacency 
matrix of the Dynkin graph and all its components are positive. However, we can 
impose the condition ( |3.4] ) in a weaker sense as an analytical continuation from positive 
integer values. 

As the observables (|3.5| ) are associated with the M-matrices, we can integrate out the 
C-matrices in the partition function. After the integration the partition function ( |3.2| ). 
can be written, up to a factor equal to the volume of the symmetry group ®x^xSU{Nj:), 
as an integral with respect to the eigenvalues A^^^;, i = 1, A^^, of the hermitian matrices 

/ nn---^--^-*-^- (3.0. 

The explicit dependence on the second coupling k can be eliminated by a linear change of 
variables 

^ iM, + 5, g,^ K^g,. (3.7) 



After that the partition function (^^) takes the form 

where the cubic potential depends only on the temperature T, 

Vx{z) ^ SxV{z), Viz) = + ^z' - ^z\ (3.9) 

The coupling k, the lattice cosmological constant, reappears in the definition of the planar 
limit, originally given by (|3.4D , which now becomes 

gsNx ^ ^Sx, gs ^ 0. (3.10) 

The eigenvalue integral ( |3.8| ) is not convergent. To make it convergent, we will drop the 
absolute value in the denominator and will understand the integrals over the eigenvalues as 
contour integrals in the complex plane. The new integral will have the same perturbative 
expansion as (|378|). We choose the integration contour C so that it starts at Xi^x = — oo, 
passes near the origin and then goes to infinity along a direction in which the potential 
grows (Fig. 5a). This is a standard prescription used to define rigorously the matrix 



integrals that appear in 2D quantum gravity, see e.g. |19]. To avoid the poles of the 
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integrand at X^^i = —Xx',j we also require that the contour C does not intersect the 
contour C obtained from C by a reflection z —* —z. 





















c 




c \ 












a 




Fig. 5a. Integration contours in the complex plane of the eigenvalues for cubic potential 
Fig. 5b. Integration contours relevant for the quasiclassical limit 

There are several choices for the second branch of the contour C, which however lead to 
the same perturbative expansion. 

3.2. The matrix chain as a chiral conformal field theory 

It is known that a large class of matrix models can be reformulated as conformal fleld 



theories |2g] flJ] p| [|3[ |T|] |2J . The correspondence MM^CFT can be estabhshed if the 
measure in the eigenvalue integral can be produced by screening charge operators. We will 
show that this is the case for the integral ( |3.8| ) where the integration goes in the complex 
plane along the contour C. 

We associate with each point x G X a chiral boson ip^iz) two-point correlator 



{0\iP:^{z)iP:^'{z')\0) = S^^^f \0g{z - z'), 



(3.11) 



where the right /left vacuum is annihilated by the negative/positive power part of the 
Laurent expansion of the flelds dipx{z) in z. We introduce two more bosonic flelds, 



Vx{z) = fx{z) - (fx+i{-z), 
and the conjugated fleld deflned by 



(3.12) 



z — z 



(3.13) 



The fleld ^x{z) is related to ipx{z) and T^iz) by 



(3.14) 
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2$,(;^) - J] = r,(z), (3.15) 

x' 

where A^j./ is the adjacency matrix of X. 

It is straightforward to check that the integrand of ( |3.(j| ) can be written as Fock space 
expectation value of a product of vertex operators 

E^{z) =: e-^-(^^ : : e^^+i^^^ : . (3.16) 

Let {N\ be the charged vacuum state defined by the asymptotics 

d^^{z) = — ^ + ---, z^oo. (3.17) 
Then using the OPE of the vertex operators one finds 



(3.18) 



To generate the external potential it suffices to replace the left vacuum by a coherent state 

{B\ = (iV|exp(Vif,), H^ = - i -^V^{z)d^^{z) (3.19) 
^ 9s J 27rz 



where the contour Coo encircles both C and C (Fig. 4a). This is easily verified by applying 
the operator product expansion ( |3.13| ) to commute e^^ to the right. The result is a factor 
e~9^^^ per eigenvalue. 



To write the Fock space representation of the eigenvalue integral ( |3.8[ ) we define the 
screening operators 



= J^dzE^iz). (3.20) 
Then we can write the eigenvalue integral in the form of a scalar product 

Z = {B\n), (3.21) 

where {B\ is defined by ( p.l9| ) and 

|0> =exp(^Q,)|0>. (3.22) 

a; ex 



• Correspondence between the observables in CFT and the matrix model 
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The expectation value ( in the matrix model becomes the vacuum expectation 
value in terms of Fock space. To any observable O^m matrix model one can associate 

an operator O such that 

(B\o\n) 

(^mmU = ^^^-((^)- (3.23) 
In particular, the resolvent of the matrix M^; is represented in the CFT as 

W^x(^)=^7.tr(^— ^ -<7sa$L+l(^), (3.24) 

where d^^^^z) is the negative part of the Laurent series of 9$(z), which annihilates the 
right vacuum, d^^^^\0) = 0. 

In other words, the resolvent is equal to (minus) the singular at ^ = part of Laurent 
series of d^{z), which annihilates the right vacuum. 

The other important observable is the FZZ-brane, which measures the effect of sup- 
pressing the integral with respect to one of the eigenvalues of the matrix at assigning 
to it a complex value z. This matrix observable corresponds to the vertex operator ( |3.16|) , 

e-i^==(^)det(z-M^)2 JJdet(^ + M^O"^"^' ^ E^{z) e'^^^'^ : . (3.25) 

x' 

In the limit ^ oo we can replace ( : e~^'^^^^ : ) = '"^'^^^ ^. Thus the expectation value 
of the field ( p.l5|) gives the effective potential of one eigenvalue taking value z in the mean 
field produced by the other eigenvalues, 

-K(2:) -2trlog(2-M,) + VA,,4rlog(z + M,0 ^ T,{z). (3.26) 

• Virasoro constraints 

In the matrix model formulation, the n-loop amplitudes satisfy a set of loop equations 
that follow from the translational invariance of the matrix integration measure. The loop 
equations have been first derived from the combinatorics of planar graphs |p. They have 
the form 

W^zf + ^f -^WAO\J2Axx'W^{-0-J2^x{0]) =0. (3.27) 
2txiJc z-Q ^ J /mm 

The integration contour C_ separates the eigenvalue contour C from its reflection image C 
(Fig. 4). It is also assumed that the point z is outside the integration contour. 

In the Fock space representation, the loop equations are equivalent to the following 
set of operator identities satisfied by the right vacuum 

i]^£^Mli |0) = o + ... (xgZ+). (3.28) 



13 



The omitted terms vanish after multiplying with the left vacuum {B\. In other words for 
each X the singular at 2 = part of the mode expansion of the energy-momentum tensor 

T,{z) = \ ■ [dT^{z)f : (3.29) 

annihilates the state |0). 

Proof: The operator {dTj^Y ^■cts to the right vacuum as 



z + C 

Jc c", ^ ^ s ~7 Jc 



(3.30) 



(^ + 



The first term disappears after multiplying with the left vacuum since the potential is 
chosen so that it grows to + infinity at the endpoints of C. The second term is regular in 
the domain inside C_ and therefore does not contribute to the integral ( |3.28| ). 

• Analytical properties of the collective field 

Arguing as in 0, we expect that in the quasiclassical limit Qs ^ the eigenvalues 
of the matrix M^; condense on some interval [—A, —M] on the negative real axis. The 
endpoints of the interval depend on the couplings T and k. Up to exponentially small 
terms the contour C can be restricted to the negative real axis, and the contour C - to the 
positive real axis. Then the contour C_ goes from —00 to below the real axis and then 
back to —00 above the real axis (Fig. 5b). 

The classical value of the resolvent 

W^iz)=gs(triz-M,)-') (3.31) 



is a meromorphic function with a branched cut along the eigenvalue interval and asymp- 
totics at infinity as 

W^{z)^gs-^ = S^-. (3.32) 



z z 

2 



The classical Virasoro constraints ( 3.28 ) imply that (dTj;) is analytic in the vicinity of 



the interval [—A, —M]. Since dTx is discontinuous across the cut, this means that 

gsdVliz) = dV,{z) - 2W^{z) -J2A..>W^,{-z) (3.33) 

x' 

satisfies the boundary condition 

dVliz + iO) + dV^iz ~ iO) = 0, ze[-A,-M]. (3.34) 
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The boundary condition (|3.34|) , together with ( p.32|) and (|3.33| ), determines uniquely the 
meromorphic function 



Fig. 6: Riemann surface of the classical solution VF'^(z). In the first (physical) sheet 
there is only one cut [—A, — M] while all other sheets contain two symmetric cuts. 

With the particular potential ( |3.9|) , the classical solution must be of the form 

K(z) = ^^'i^), ri{z) = ^T^{z), W^{z) = ^W^{z) (3.35) 
oi oi oi 

and ( |3.34| ) becomes a condition for the single function W'^{z): 

W''{z + iO) + W%z-iO) -2cosnpoW''{-z) ~dV{z) =0, ze[-A,-M]. (3.36) 

This is the same equation as in the 0(n) matrix model with n = 2cos(7rpo) [|25|. The 
solution of ( pj.34| ) for a general polynomial potential can be expressed in terms of Jacobi 



theta functions . Its explicit form in the case of the cubic potential ( |379| ) can be found 
in 



2% . The Riemann surface of the function W^{z) is sketched in Fig. 6. 
It is consistent to assume that the quantum field has the same analytic properties as 
the classical solution and define its mode expansion at infinity in terms of a complete set 
of functions globally defined on the Riemann surface. This is the basic assumption for 
the iterative procedure for calculation of the n-loop functions perturbatively in Qs, known 
as method of moments, Our method can be considered as an operator equivalent 



of the method of moments. We thus assume that the boundary condition (|3.34|) actually 
holds at operator level, 

d,[^^{z + iO) + $,(z -iO)-J2 ^x,x'$x'(-2)] = [r,(z + iO) + r^{z - iO)] 

(3.37) 

= 0, 2G[-oo, -M], 

and should be understood as a condition on the n-point correlation functions of the col- 
lective field $. 



2 



Note that the boundary condition ( 3.34| ) has a natural interpretation in the matrix model. 



It means that the classical distribution of the eigenvalues is such that the effective potential 
T%(z) is constant on the eigenvalue interval. Since the derivative of the effective potential has a 
discontinuity across the cut, one has to take the half-sum of its values on both sides. 
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3.3. Continuum limit, phase diagram and critical points 

In the ensemble of dynamical triangulations the critical thermodynamics is controlled 
both by the temperatore T and the bare cosmological constant k. The continuum limitl is 
achieved, for given T, at some critical value Kc{T). The singularity along the line k = Kc{T) 
is due to the contribution of surfaces of diverging area. In the continuum limit, i.e., near 
the critical line k = Kc{T), the lattice can be approximated by a continuous worldsheet 
and the theory is described by perturbed Liouville gravity. 

The behavior of the theory in the continuum limit depends on the value of T. As in 
the case of a flat lattice, there are three possible types of critical behavior, characterized 
by massive, dilute and dense loops |^[^. The critical line k = k,c{T) consists of two 
branches: 

( k[{T), if T> T*; 
[Ki^iT), i(T<T*. 

The first branch k = describes high-temperature phase where the area of the graph 

diverges, while the loops remain finite. Its equation is given by dn/dM = 0. Along 
this branch the critical behavior of the partition function is that of pure gravity with 
Cmatter = 0. The secoud branch k = k,ii{T) describes the low-temperature phase of densely 
packed loops. It is determined by M = 0. Here the critical behavior is that of a 2D gravity 
with matter central charge given by (|2.14|) . The two branches meet at the critical point 
K — K* ,T — T*, where both the area of the graph and the length of the loops diverge. The 
double scaling limit at this point describes the phase of dilute loops with matter central 
charge (|2.13|) . 



The vicinity of the critical point T* , k* is parametrized by the renormalized cosmo- 
logical constant fx and the temperature coupling constant t, defined as 

~ - K, t^T* -T. (3.38) 

We will consider the more general case of finite t, where the matter field is not conformal 
invariant. Then the susceptibility u = —4pod'^Z, where Z is the partition function on the 
sphere without boundaries, satisfies the transcendental equation P7[| B 

^=i^^ + 23fei«^- (3-39) 
This equation describes the flow between the dilute {t = 0) and the dense {t — > oo) phases. 
The dimension of the coupling t ~ matches the gravitational dimension Si^s = 1 — po, 



obtained from Ai^3 by the KPZ scaling relation A = ^^^^^^^ pO[| . This is consistent with 
the conjecture that for flnite t/iJ,P° the theory is described by a perturbation of the critical 
point t = by the Liouville-dressed thermal operator $1,3. Eq. ( |3.39| ) is discussed in the 
context of perturbed Liouville gravity by Al. Zamolodchikov ||31||. 



^ We do not call it scaling limit because the theory is scale invariant only at the critical points. 
^ The normalization factor 4po in the definition of the susceptibility is fixed by the comparison 
with the 3-point function, see below. 
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3.4- The loop operator and the disk partition function in the continuum limit 

The critical line k^J (T) is determined by the condition that the right edge of the eigenvalue 



distribution (after the shift (|3.7|) ) hits the origin while the left end stays at a finite distance 
A ~ T. Therefore in the continuum limit M <^T. If we introduce the renormalized length 

I = ^L, (3.40) 
then the loop operator trM^ can be approximated, after the shift (p77|), by the exponential 

tr[(M, + iT)^] ^ e^^'^^°sT^^^g£M.)_ ^3 41) 

Therefore up to trivial multiplicative factor the operator creating a boundary of length i 
and height x is 

W^.(£)=^7str(e^^^). (3.42) 

The Laplace transform 

= / d£e-''W^{i) = gM ^) (3-43) 

is the operator creating a boundary with a marked point with height x and (in general com- 
plex) boundary cosmological constant /is = z. The n-loop amplitude ( |3.5| ) as a function 
of the boundary heights Xi and the boundary cosmological constants Zi is given by 

Z(xi, zi; . . . ; Zn) = Jim^^7s"" ^ f] tr log {z, - M,J ^ . (3.44) 

4. The collective field theory as a CFT on a Riemann surface 
4.1. The classical solution in the scaling limit 

In the following sections we will show how Virasoro conditions ( p.28|) can be solved directly 
in the continuum limit A — 00. The solution is completely determined by the classical 
value of the resolvent, whose explicit form we give below, and the target space X of the 
height model. From now on we will concentrate on the SRSOS model, 

X = Z_|_, S^^^ = sin( Tcpx), = sin(7rpo3:^), V'p(^) = " / ^ 1 . 

sm(7rpoa;) 

We will give an operator solution of the Virasoro conditions ( |3.28| ). The n-loop amplitudes 
will be expressed as Fock space expectation values as in (|3.23|) , but the left and right 
vacuum stated {B\ and |0) will be given a new realization in terms of a chiral bosonic field 
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on this Riemann surface. The basic idea of this approach is that the coUective field can 
be constructed perturbatively in the string couphng Qs as a bosonic field defined on the 
Riemann surface of the classical solution [ |TB| , P^ . The operator solution yields a diagram 
technique similar to the one derived for the ADE matrix models . 

In the continuum limit the classical solution is written in terms of the uniformization 
parameter r by 

z = M coshr, 

1 f M'^+P" M'^-Po \ 
d^^''\^ = - — [— cosh(H-po)r -t-^ cosh(l-po)^ , (A-\) 

1 MPo 

d^,^''\z = coshpoT- 

^9s Po 

The expressions for the two derivatives are compatible if they lead to the same expression 
for the second derivative. From here we find, using the formulas 

= d^M {Om - MtLhr ^r) , = M^dr 

a transcendental equation for the modulus M: 

2ij = + j^tM^-^P'>. (4.2) 



It can be shown that the modular parameter M is related to the string susceptibility 
u by 

m = M2po, (4.3) 

which relation implies the equation of state ( |3.39| ). 

From the solution ([4.1|) we obtain for the classical effective potential 



d,T''{r±in)=±id,ip''{r), (4.4) 

where 



a^(^^(r) : = 2sin(7rxa^)a^$^(r) 

smh(l + po)t + smh(l - po)r 



(4.5) 



This function has two cuts, [— oo, — M] and [M, oo]. The two signs in (|4.4| ) correspond to 
the value of dV^ above and below the left cut. 

The field 99'^ has vanishing imaginary part along the interval z > M, which implies 
that dT'^^{z) has vanishing real part along its left cut z < —M. Therefore it is expanded at 
the branch point 2 = — M as a series of the half-integer powers of 2 + M. In the following 
we will need the explicit formula for this expansion. 



d,V%{z) = S^y , 1 + V7) • (4-6) 

^ M^g, ^^^(2n-l)!! \ m) ^ ' 



18 



The coefficients Hn, n = 1,2, wliidi we call moments of the classical solution, are given 
by I 

= 2 (m2+poF^_i(1 +po) + tM2-P0F„_i(l -po)) , (4.7) 

where 

FM ■■= (4.8) 
n! 

and {y)n is the Pohgammer symbol. In the particular case po — (pure gravity) the series 
1) truncates to the first two terms. 



4-2. The collective field as an operator field on a Riemann surface 

In the parametrization z = M cosh r the operator boundary condition ( |3.37| ) takes the form 
of a difference equation with respect to r. Considered as a function of the uniformization 
variable r, the operator field ^xij) = ^x{z)\z=m coshr must 

i) be entire function of r, 

ii) satisfy the orbifold condition 

$x(r) = $.(-r), (4.9) 

iii) satisfy the discrete Laplace equation in the {x, r} space: 

x' 

The last condition can be conveniently written as a difference equation 

(cos Tr^r - cosh a^) $^(r) = 0, (4.11) 

if the restriction to positive heights is imposed by requiring antisymmetry = —^-x- 
The condition ( [4.11|) means that the components {^x{z)}x^z+ and {Tx{z)}x^i,^ of 



the collective field can be obtained by analytical continuation from a single holomorphic 



^ The expansion coefficients can be obtained from the standard representation of the hyper- 
geometric function (we use that \/ ~ '^sinh ^-^): 



ir\\2 "T i/' 2 y' 2' 2M I ~ sfi (2n-l)!! \^ ^ M > 



i sinh(g(T — itt)) / z + M 

2g — y 2M 

n>l 
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field $(t) := $i(r) . Indeed, ( |4.1C1| ) is equivalent to the defining functional relation of the 
Chebyshev polynomials of second kind 

$x(r) = U^-i{cosndr) $i(r) 

sin nxd-r ^ , , (4-12) 

= -■ TT- ^{r). 

sm -kOt 

For the field '^xij^ we have, 

rjr ± m) = 2$,(r ± in) - $.+i(r) - $:,-i(r) 

= ±2zsin7r9T- $x(t) = ±2ismnxdr $(r) 

= TT- ^{r) 

sm Tro-r 



(4.13) 



for real r. Here we introduced for convenience a second field f related to $ by 

ipir) = 2sin(7ra^) $(r). (4.14) 
Note that $ is even, $(r) = $(— r), while is odd, v^(r) = — (/?(— r). 
^.5. Mode expansion and two-point correlator 

We write the operator field $ as a sum of positive and negative frequency parts, 

$ = $+ + $-, (4.15) 

so that $~ as a function of z is analytic on the Riemann surface with the point z = oo 
removed, and is analytic in the vicinity of the point ^ = oo. The two pieces are given, 
as functions of r, by the spectral integrals 

oo 

nrvaVi y/T ff) 

1/ 5 



^ (''")= / 17 coshj^r 

.1 (4.16) 
$+(r) = - / ^ e"^^ $+, Re(r) > 0. 

Similarly, for the quantum field ^{t) we find, using ( [4.14| ), the mode expansion 



oo 

2 sin -KU 



(^-(r)= / du sinz/r$;, 

.oo ' (4.17) 

The operator amplitudes are assumed to satisfy the canonical commutation relations 

[$+,$;,] =z/5(z/-z.') (4.18) 
20 



and the left and right vacuum states are defined by 

(Ooo|$-=0, $+|0oo) = 0. (4.19) 
With this definition of the left and right vacuum states the field has zero expectation 



value. The expectation value (|4.1| ) can be generated by replacing the left vacuum by an 
appropriate coherent state. The two-point correlator is given, for Rer > Rer', byi 



(0^|$(r)$(r')|0^ ) = ilog(r2-r'2), 

(0^ \'^{r)(f{r')\0^ ) = 2sin7rar sinTT^r' log(r^ - r 



(4.20) 



where we neglected an infinite constant term. From here we evaluate the two-point func- 
tions of the original field and ^x'- 



:*.M*..(r')) = (o„ H.(r)t(r')|0^ ) 

sm^TrOr) sm^TrOr') 

x—l x' — l 

= J2 Yl (0^ |$(r + Z7rm)$(r' + W)|0^ ) 



(4.21) 



m=l — x m' = l — x' 



sm{ndr) sm{ndr') (4.22) 
= -4sin(7ra;a^) sin(7ra;'a^) (0^ |$(r)$(r') |0^ ). 



The two-point correlation functions (|4.21| ) and ( [4.22| ) are diagonalized in the momen- 
tum space, 



$(r, p) = sin yrpa; $a;(r) , if{T,p) = ±- sin Tcpx T x{t ± 



ITX . 



(4.23) 



xEi 



The two fields are related by ( f4.14| ), 

(^(r, p) = 2 siuTT^r $(r, p). 



(4.24) 



Performing the sum in x in ( [4.23|) we obtain the mode expansion of $(t, p) and (/?(r, p) , 
which is a restriction of the spectral integral ( [1.16{ ) to a discrete sum over the momenta 
u = ±p(mod2). For ip{r,p) we find 



= E {%+2n\ sinhdp + 2n\r) + $+ 



P+2n\ 



\p+2n\T 



(4.25) 



^ The normalization of the field was fixed by the condition that short distance behavior at 
T = r' is compatible with ( |3.11 ) and ( 3.14| ). 
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The two-point function of Lp{p, r) is, for Re(r — r') > 0, 



0„ !»>(?, rMp',T')|0„ ) 



-\p+2n\{T-T') _ g-|p+2n|(T+r') 

\p + 2n\ 



P — P 



(4.26) 



Here and below 5^^^^ denotes the periodic (5-function 

5(2) (p) :=J2S{p + 2n). 



(4.27) 



The time-ordered correlator 



d^^\p-p') G{T,T'-p)=e[Re{T-T')] (0^|vp(p,r)(^(p',r')|0^ ) 

+ ^[Re(r'-r)] ( 0^ |^(p', r')^(p, r)|0^ ) 



(4.28) 



is further diagonalized as 



/>oo 

/ drdr' sirv^Er) sin(£;V') ^(r, r'; p) = 2tt 5{E - E') G{E, p), (4.29) 
Jo 



G{E,p) = J2 



1 



TT sinh Tri? 



1 



ne2 



+ (p-K2n)^ 



i? 2 cosh TT-E — 2 cos np 



(4.30) 



Finally, using the relations ( [4.23| ) and ( [4.24|) we find the spectral integral for the the two- 



point correlator (|4.21 



($.,(ri)$.,(r2)) 
A{E,p) 



71 



dE / dp cosEti cosEt2 siuTrpxi sm7Tpx2 A{E,p) 
Jo 

n 1 



(4.31) 



E sinh 71 E 2 cosh 7rE — 2 cos 7vp 



This is the well known expression for the annulus amplitude in momentum space 
3^ , p3| ,|9[] . The collective field ^x{t) creates Cardy type boundary state |t, x). The in- 



termediate states in the spectral integral are Ishibashi states are characterized by two 
quantum numbers, the Liouville energy E and matter momentum p. 

The spectral decomposition ( [4.31|) is valid also for ADE string theories, with the 
integral over momenta replaced by a discrete sum over Coxeter exponents. The annulus 
amplitude for the A-series, or RSOS models coupled to gravity, was derived using the world 
sheet CFT approach in [0, see also for an expression analogous to ( [4.21|) . 
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5. Operator solution of Virasoro constraints on the Riemann surface 

The idea of the operator solution of Virasoro constraints (|3.28D is that it is possible to 



introduce another Fock space, associated with the the mode expansion of the fields Tj:{z) 
near the branch point z = —M . Then the problem reduces in the scaling limit to the 



already solved problem of perturbed Cmatter = — 2 gravity where the operator solution 
of Virasoro constraints has been given in terms of a twisted boson. In this section we first 
define the mode expansion of the field F^; in the half-integer powers of z + M and the 
associated bare Fock vacua, (0, I and 10, ) Then we write down the solution of Virasoro 

' \ tw I I tw / 

constraints in terms of the corresponding oscillator modes. This will lead to the operator 
replacing the right vacuum Finally, we will construct the operator replacing the 

coherent state {B\. As a result we will obtain a diagram technique similar to the one 



obtained in 141 for the ADE models. 



5.1. Mode expansion at the branch points z = —M 

In the following we will use also the rescaled string interaction constant and moments 
defined by 

Qs = — , Hn = — • (5-1) 

Hi III 

For each x, the operator field rj;(z) behaves near the branch point z = —M as a twisted 



2 

boson. Its mode expansion is given by a series of half-integer powers of z+M = 2M cosh ^ 
now positive and negative. 



2 ' 



n>0 



(2n- 1)!! 



(5.2) 



n>0 



In this expansion we dropped all terms in the expansion of the classical value ( [4.6| ) except 
the first one, proportional to The rest of the classical background will be later 

reintroduced as a perturbation. 

The operators an,x and aj^ ^ and the twisted left and right Fock vacua satisfy 

(0,JaJ,,, = 0, a^,,|0,J =0 (n > 0, x G Z+). (5.4) 



^ If we had dropped the whole classical field, Virasoro constraints would become singular and 
would have no solution 
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The vacuum state |0^^) can be thought of as the direct product of an infinite number of 
twist operators associated with the branch points of the Riemann surface [pT] , ^ . The 
two-point function of the field rj:{z) in the twisted vacuum is 

, V cosh 5 — cosh ^ 

(0,jr.(z)r,,(z')|0,J =5,,x' In— -Jr. (5.5) 

cosh ^ + cosh Y 

The creation operators ^ can be expressed by the generating function 



oo 



where the contour C_ encircles the cut [—00, —M] of rj:{z). Due to the vanishing of the 
real part of F^; along this cut the contour reduces to a closed circle around —M. 

We will actually need the expansion of the field (f{T,p), which is written in terms of 
the Fourier transformed creation and annihilation operators 

flfcb) = X] sin{7rpx) ak,x 

(5-7) 

= Yl sin(7rpa;) a^, 

where the canonical commutation relations have the form 

[an{p),al{p')]=6^^^,S^^\p-p'). (5.8) 

Then the twisted propagator in momentum space is given by replacing Sx,x' 6^'^\p—p') 
in (|3p. 

The formula ( |5.6|) can be used to invert the expansion ( ^.2|) for the creating operators, 
see Appendix B. We have in momentum space the operator identity 



at 



ip)-5k,i 9:'5^^Hp-po) = Fk{dr)d^ip\r,p) (5.9) 



T = 



where 99^ is given by the Fourier transform of the first sum in ( |5.2| ) and the differential 
operator Fk{dr) is defined by (^?8|), with g ^ dr, 



We shall also use a weaker version of ( p.9| ) also derived from ( |5.6| ) in which the negative 
mode part in ( [4.15| ) appears in the r.h.s., 

{a\{p)-5k,i9:H^''\p-PQ))\Qoo) = Fk{d^)d-rip-{T,p) |0oo). (5.10) 

T = 
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5.2. Operator solution of Virasoro constraints on the Riemann surface 

The right physical vacuum state wiU be constructed from the modes of the Z2-twisted 
boson. The conformal invariant physical vacuum \Vt) must be of the form 

1^^)= n (5.11) 

where the operator Q^; represents a formal series expansion in terms of the creation opera- 
tors aj^ 3.. The coefficients of the series are fixed by Virasoro constraints ( p.28|) , where the 
stress-energy tensor ( p.29| ) is expressed in terms of the oscillator modes associated with 



the expansion (p.2|), 

T^{z) = is[ar.(^)]2s = ^^L^,.(l-f^)— 2, (5.12) 

n 

where □ S means the normal ordering with respect to the mode expansion ( ^.21) . We write 
down explicitly the first two of the generators L^, n > —1, 



Lo,x = ^ ^^(2m + l)(aJ„^^-(5^,i|^)a^,^j +^ 



(5.13) 



We are looking for a perturbative in solution of the conditions of conformal invari- 

ance 

L,,, |O) = 0, n>-l, 

of the right physical vacuum state ( |5.11|) . The problem is identical (for each x) to the one 
of solving Virasoro constraints in pure gravity. The spectrum of operators in this theory is 
described by the flows of the KdV hierarchy [p6| , |37[| . An important property of the theory 
is that there is a conserved charge (ghost number) with a distinct background charge for 
each genus, so a specific correlation function can be non-zero for at most one genus. The 
explicit form of the operator is 



(5.14) 



where (r^^ ■ ■ -tj.^ )^ is the genus g correlation function in the topological (cmatter = —2) 
gravity. 
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Note: The generating function of these correlation functions can be written as 



Sx 

9s 



1 
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Sx 

9s 



exp Yl 

\n,g>0 

2-2g-n 



Sx 

9s 



2-2g-n 



E 



(5.15) 



/ci + ...+A;„=n+3g — 3 

QU 



logZx{t)i^__ 



(5.16) 



(^tw • • • ^hji ^~^X | ^tw ) | conn 

The Virasoro generators are reahzed as differential operators, i.e., clI, ^ ^ tk , dk,. 
The formulas in |3^] differ from ( ^.13|) and (|5.15|) , (|5.16D by the rescaling Sx/gs 1- 

We will restrict ourselves to the genus zero correlation functions and consider only the 
piece with = in the sum in ( ^.14|) , for which [^ , |38|| 



I \ (mi + ...+m^)! 

(Tmi---rm„)o = i i , mi + ... + = n - 3. 

m\\...mn\ 

In momentum space we have, neglecting the overall constant, 

1 

|0)=exp (^I"^^^"^ y dvi...dv^Y Vfc,...fc„ (pi, 4^ (pi)....4^(p. 



(5.17) 



n>3 



A;i + ...+fc„=n— 3 



with 



Vfc,...fc„ (pi,...,p^ 



(/Cl + ... + /Cn) 



|0.J 
(5.18) 

(5.19) 
(5.20) 



iVpi,...,P„ = XI '^^^ ^Pi(3^)-V'p„(a:^)- 

5. 3. Classical background and two-point correlator in terms of twisted bosons 

The gaussian field 9rj;(r) is completely characterized by its vacuum expectation value 
and its two-point correlator ( [1.22|) . In terms of the mode expansion ( [4.17|) , the left physical 



vacuum {B\ in (|3.21| ) can be expressed as the coherent state associated with the vacuum 
(Ooo| that generates a vacuum expectation value. Now we want to express it in terms of 
the mode expansion (|5.2|) and construct {B\ as a perturbation of the left twisted vacuum 
(Oj^|. In this case the perturbation is not just a coherent state but rather Bogolyubov 
transformation, which changes also the two-point correlator. 



(S| = (0,Jexp(l Yl E 



k,k'>0 x,x'e'L+ 



Sx \^ \^ ^ ri 

O'k'x' ~ / , / , l^k^xCtk,x 

Qs 

k>2xe1+ 



(5.21) 
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In momentum space with 



S'^'^\p-p') Dkn{p) = ^ sm{npx) sm{'Kp'x') D^'^ (5.22) 
( ^■21|) is rewritten as 

1 

{B\ = {0,Jexp(^- dp Dk,k'{p) ak{p)ak'{p) - — Yflkak{po)y (5-23) 

k,k'>0 _^ k>2 

The left vacuum defined in this way has the properties 

{B\ {al^-Sn,i^)\0,J = -ftn^, (5.24) 

9s 9s 9s 

The first of these identities is equivalent to 

{B\T^{z)\0,J = {0^\T,{z)\0^)+T:{z) = TI{z) (5.26) 

The unknown kernel D^'^, in the second identity ( ^.25| ) is determined requiring that i 



(s|r,(^)r,.(;^')|o.J = {o^\r^{z)r^,iz')\o^) + rUz)K,{z'). (5.27) 

Prom the expansion (|5.2D we obtain 



{B\TUz)rl{z')\o,j = (Ooo|r.(z)r.,(/)|0oo) - (o,jr.(^)r.K^')|o.J + r^(^)r^.(^') 

(5.28) 

where the twisted 2-point function is given in ( ^.5|) . 

To calculate the kernel Dnn'{p) in momentum space, we apply the operator identity 
( ^.9|) to the Fourier image of ( |5.28| ) and compare with the Fourier image of ( |5.25|) . As a 



result we obtain an expression for Dkk' {p) in terms of the difference of the 2-point functions 
in ( p.28| ). For the Fourier transform ( |5.22| ) the relation reads 



Dkk'ip) = drdr'Fk{dr)Fk'{dr')D{r,r',p)\r=0=r' , (5.29) 



^ This identity as well as (|5.26| ) is consistent with the quahtative identification of the left 
vacuum {B\ with (Ooole'^*^ , where the positive mode ^t=g generates the classical part $'^(r) of 
$^(r). On the other hand we can identify the right vacua |Otw) and |0oo)- 
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where 



It is useful to note that the second term in ( ^.301) is equal to the first term at p 



(5.30) 



(O.J<^(r)vp(r')|0,J = G(r,r',p=i) 



sinh ^ — sinh ^ 

In 4- 

sinh ^ + sinh ^ 



(5.31) 



The derivative of the first term in ( |5.3UD computed from ( |4.26D takes, up to the (5-function 
factor, the simple form 



dr{0^\^{r,pMr',p)\0^) 



coshdpl - l)(r + r') cosh(|p| - l)(r - r') 



sinh(r + r') 



sinh(r — r') 



(5.32) 



The difference in ( p. 301) is smooth everywhere and can be expanded in Taylor series at the 
point r = r' = 0: 



D{r, r'-p) = dr {\p\ - i)(|p| - I)tt' [l + ^{p^ - 2\p\ - 1){t^ + r 
+ Up'-Mp\' + y + ^IPI + i) (^rV- + ^(r^ + r-)) + ■ ■ 



(5.33) 



e.g., 



2^ 



i^oo(p) = (N-i)(H-i), 

DM = -\i\v\-\m-im + mp\ 
Diiip) = i{\p\ - km - l)i\p\ + i)i\p\ - Dip" - 2p - !), 

5.4- Diagram technique for the correlation functions 



(5.34) 



etc. 



Now we are in a position to compute of the n-loop amplitudes. The scaling limit of the 
amplitudes ( |3.44| ) is described as in ( p.24|) by the positive mode part of the field $(^). 
Thus the n-point function for n > 3 will be defined as 



($(ri,pi)...$(r^,p^)) 



genus zero 



^72-"(S|$+(ri,pi)...$+(r„,pOI^^)- 



(5.35) 



The field $''~(r, p) is the Fourier transform of the loop creation operator which is assumed 
to annihilate the twisted vacuum $+(r, p)|Oj^) = 0. 

In order to evaluate the n-point amplitudes (|5.35|) we will also need the commutation 
relation between the loop operator $"*"(!), r) and the creation operators aj^. This can be 
done using the linear relation ( ^.10|) between aj. and ip~{p,r) and the expression for the 
commutator of <f>~^ and which is equal to the two-point correlator ( |4.26|) . Using the 
expression (|5.32| ) for the r-derivative of the propagator, we find 



Ooo|v'(^, P)ai{p)\Oc 



{k-dr)k{k+dr)k cosh(l - \p\)r 



k\ 



sinh r 



(5.36) 
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Then from the relation ( |4.24|) we get 



We would like to write the leg factors directly for the local operators obtained by 
shrinking the boundaries to punctures. For this purpose we expand the loop operator 
as a power series in z. The allowed powers are of the form z~'^ with u = \p + 2n|, n G Z. 
The corresponding amplitudes create local operators on the world sheet with gravitational 
scaling dimensions ^{i' — po)- Therefore the leg factors for the local operators are given by 
the coefficients in the expansion of the leg factors for loops, 

Fk{p,r)= Yl Fk{u){z/M)-\ (5.38) 

u=\p+2n\,nel' 

The coefficient corresponding to = p is obtained by replacing the last factor on the 
r.h.s. of (|5.37| ) with and substituting the derivative 9^ by u'^ = p^. Thus the n- 



tachyon correlation function ...,Pn) is obtained from the n-loop amplitude ( p.35|) by 

replacing 

These are the leg factors for the order operators with \p\ < 1. 

Now we are able to evaluate the n-point function G{j)i^ ■■■,Pn) of order operators by 
performing the necessary number of Wick contractions in ( p.35|) . This amounts to a sum 
of Feynman diagrams composed by vertices, propagators, tadpoles and leg factors. The 



leg factors come from the commutators ( |5.37|) , while the tadpoles are associated with the 
action of on {B\ as in ( |5.24| ). The tadpoles are proportional to the normalized moments 
An = /^n//^i and depend on t and /i through the dimensionless coupling 

i=^7^- (5-40) 
The first moment /ii, eqn. (f4.7|), then reads 

/xi = 2M2+po(1 + £). (5.41) 



Since the string coupling constant enters only through the ratio Qg = Qs/^i, the terms 
retained in the perturbative expansion of an n-point function are only those with an overall 
factor The vertices can be attached either to tadpoles or to the propagators or to 

leg factors. We summarize the Feynman rules in Fig. 7. 
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Leg factor 



Propagator 



Tadpole 




^1 P2 



Fig. 7 : Feynman rules for the correlation functions 



Feynman rules: 



External line factors (leg factors) Fk{p), where 



Fkip) 



Propagator Dk,k'{p)- 



Tadpole Bk{p): 



Dooip) =i\p\-im-i) = -Fiii-\p\). 

Doi{p)=-F2il-\p\)=Dio{p), 
D^,{p)=-2F2il-\p\)-2Fs{l-\p\), . 



Bo{p)=Bi{p) = 0, 

Bn{p) = -S'-'^HP-Po) An; 

Fr,-l{l+Po) + t Fn-l(l - PO) 



1+t 



n>2 



Vertices Vk^^...,k„{Pi^ —^Pr 
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• An overall factor 

gl^-'^Y[MP' = g^-^ ;Ui2-"MPi+-+P". (5.46) 

i 



We derived the Feynman rules for the case of the SRSOS model, where the momenta 
belong to the interval < p < 1. One can derive a similar set of Feynman rules for the 
n-point functions in the SOS model, where the the momentum interval is — 1 < p < 1. In 
both cases the multiplicities A^pi,...,p„ are given by 

Np,...,p. = Yl ■ ■ ■ (5-47) 

where si^^ is defined by (p.l2| ). 

The multiplicities ( |5.47|) are periodic in p ^ p + 2 and, in the case of the SRSOS 
model, also antisymmetric in p —p. The propagator Dk^^'ip) is the same for both 
models. It has the symmetries Dk^k/{p) = Dk,k'{—p) = D^^k'i'^ ±p)- The propagator and 
the vertices can be defined outside the interval |p| < 1 by periodicity in p — > p + 2. 



Fig. 8: The only Feynman diagram for the three-point function 

The simplest example is the 3-point function (Fig. 8). The only term of order 
in the exponent in |0) which contributes to this correlator is the one with (toToTo), i.e., 
/ci = = /c2 = /ca. The corresponding leg factors for the local operators are -Fo(Pi) = 1 
and we obtain 

j\^Pl+P2+P3 j\^Pl+P2+P3-2-po 

^(^1,^2,^3) = Vfc,,fc2,fc3bl'P2,P3) = 1 ^ tM-2P0 ^(^'l'^2,P3) • (5.48) 

The 3-point function is proportional to the 3-point fusion multiplicity. In the case p\ = 
P2 = Ps = Po it is equal to the third derivative of the partition function on the sphere. 
Indeed, 

This justifies the normalization of the susceptibility, u = —4pod'^Z = M^^". To compare 
with the results from the world sheet CFT, we should adjust the normalizations of the 
couplings /i and t, as well as the one of the order operators. 
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6. The 4-point function in the SOS and SRSOS models 

6.1. General formula for the 4-point function 

The 4-point function is given by the sum of the three Feynman diagrams shown in Fig. 
9. Each Feynman diagram stands for the sum of terms that differ by permutations of the 
external legs. 





+ 



Fig. 9: The diagrams for the 4-point function 

We will write the corresponding analytical expression in such a form that it holds both for 
the SOS and the SRSOS models. We assume that the momenta can have both signs and 
the momentum integration is done in the interval —1 < p < 1. For the SRSOS model the 
momentum space interval folds to (0, 1) because of the symmetry of the Feynman rules 
under p —p. We have 



X 



dp [Fo{pi)Fo{p2)Np^p^pDoo{p)N-pp^p^FQ{p3)Fo{p4) + perm.] 

(6.1) 

+ /^r^ Np 

1P2P3P4 \Fo{pi)Fo{p2)Fo{ps)Fi{p4) -|- perm. 



(-A2) Fo{pi)Fo{p2)Fo{ps)Fo{p4) 
Then, using (|5.42| )-( |05| ) we rewrite ( |0| ) as 

\ s=l 

-\- J dp {Np^^p^^pN—p^p^^p^ -\- Np-^^p^^pN—p^p^^p^ -f- Np-^^p^^pN—p^p^^p^) — 2) (1^*1 ~ 2) ^ ■ 

(6.2) 

This is the general formula for the 4-point function in the continuum limit. The 4-point 
function depends on the scale parameter t = t/[M{t, /i)]^^°. It can be compared with the 
CFT results at the critical points t = and t ^ 00, where the world sheet theory is de- 
scribed by Liouville gravity with matter CFT component with central charges respectively 



Pl,P2,P3,P4 
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( p.l3|) and ( p.l4 ). At the two critical points fi2 is given by a numerical constant, 



-A2 = ^ +Po + 2^.po ^ 2+pl±2po, (6.3) 
while fxi scales as a power of M{fx, t): 

//I ~M2±po. (6.4) 

Here the sign is + for t = (dilute loops) and — for t ^ 00 (dense loops). 

In the following we will write more explicit expressions using the concrete form of the 
fusion coefficients in the SOS and SRSOS models. To simplify the algebraic expressions 
we introduce an auxiliary quantity equal to the integral of a function f{p) with respect to 
the intermediate momentum p and weighted by the fusion coefficients involving the four 
external momenta: 



)-fip)-{ ■■- 



J dp f {p) {Np-^^p^pN—pp^p^ + Np-^^p^pN—p^p^ + ■^piP4P-^—pp2P3p) • (6.5) 



In particular )— 1— ( = '^Np-^p^p^p^. Then the general formula ( |6.2|) can be written as 

]\/[\Pl\+--- + \P4.\ r -1 

G'(pi,p2,P3,P4) = -2 )- |( - /i2 + 2^ (i - P^) ) + {\P\ - i) {\P\ - f) -( 



^1 

j^|pi| + ... + |p4|-2-2po 



s=l 



2 



s=l 



(6.6) 



6.2. The 4-point function in the SOS model 



We first evaluate the 4-point function for a simpler case, that of the SOS model on random 
triangulation, whose continuum limit is supposed to be described by Liouville gravity 
having a twisted boson as a matter field. In this case the wave functions are the plane waves 
( p.20| ) and momentum-space vertex ( p.20|) is a periodic delta-function, which describes 



w(l)-type fusion rules: 

n 

Np^_p„=6^^HJ2(P^-Po) + ^Po)- (6-7) 



k=i 



Since the fusion coefficients are delta- functions, the integral over the intermediate momenta 
yields Np-^p^p^p^ times a sum over the three channels. 

We would like to avoid the complications related to the compactness of the momentum 
space. This is why we assume that the external momenta are sufficiently small, so that 
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the periodic delta functions in (|6.7|) can be replaced by ordinary ones. Then we can easily 
perform the integration in (|6.5| ), which leads to the following simple expression: 

= A^pip2P3P4 {{Pi +P2- Po)"^ + permutations^ 

4 ^ (6.8) 

S = l 

(For arbitrary external momenta we would obtain a more complicated expression, which 



is periodic in ps and therefore not analytic.) Substituting this in ( |6.6| ) we get the final 
formula, in which the quadratic in the external momenta terms cancel: 

j^|pi| + ... + |p4|-2-2po 
Gipi,P2:P3.P4) = NpiP2P3P4^ 

(1+0 (6.9) 

2 + PO - |Pl +P2 - Pol - \Pl +P3 -Po\ - \Pl +P4 -Po\ 

At the critical point t = this expression coincides with the well known formula for the 
4-point correlator of gaussian matter field []TB|] . Indeed, using the map ( |2.19| ) we write the 
expression in the brackets as 

Q - |Pi + P2 - eol - |Pi + Ps - eol - \Pi + P4 - eo\ (6.10) 

with 

Ps=bps, Q = b{2+po) = l+b, eo = bpo = I - b; 6=-=^= (t = 0). (6.11) 

b b V 1 + Po 

The world sheet effective action for the point t = is that of Liouville gravity with gaussian 
matter: 

^ j d''a[{da(P)^ + idaxf + iQ<P + ieox)R\^+Ke'''^] + ghosts, (6.12) 



where ~ fx. Therefore, at the critical point t = 0, (|6.9|) reproduces the CFT result of 
p!5[] for gaussian matter field. In the dense phase, t ^ 00, the formula ( |6.9|) can be again 
interpreted as the 4-point function with gaussian matter, but with a smaller central charge 
( P.14|) . In this case we identify 

Ps = -Psb\ Q = {2-po)b' = ^+b', eo = -pob' = l-b'; b' = (t ^ 00). 

y y VI - pq 

(6.13) 

At finite t the world sheet theory does not factorize to a matter and Liouville components. 
It can be considered as a perturbation by the Liouville-dressed thermal operator $1^3 of 
the CFT describing the dilute phase.i 



^ This perturbation is expected to lead to the periodicity of the correlation functions under 
P ^ P + 2b which matches the periodicity p — > p + 2 in the microscopic theory. We thank Al. 
Zamolodchikov for a discussion on this point. 
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6. 3. Four-point function in the SRSOS model: generic momenta 

Now we will consider the case of generic momenta in the SRSOS model. We assume again 
that the moments are sufficiently small so that the compactness of the momentum space 
is not felt. For generic (but sufficiently small) momentum p > the wave functions 

, , , sinfyrpx) 

M^) = ■ t \ (6-14) 

lead to a formal expression for the vertices ( |5.45D , which does not seem to have a rigorous 
meaning, even as a distribution. It is therefore difficult to give a meaning of the general 
formula ( |6.1| ) for the 4-point function for generic momenta. We can however modify the 
definition of the observables to obtain a sensible formula for the 4-point function. Namely 
we consider the wave function for given momentum as a sum of two terms 

^,(x) = V^W(x) - ^(-)(x), )(x) = (6.15) 

The two terms, 'i/'p^'* ^ind 'i/'p \ are analogs of the tachyons of positive and negative chi- 
rality considered in the world-sheet CFT analysis m, and are related by p — > —p. With 
this replacement we pick up the analytic expression for the 4-point function valid in the 
infinitesimal vicinity of pi = . . . = p4 = po = 0. 

In the case of irrational po we can evaluate the "multiplicities" Np^^,,,^p^ by formally 
expanding the denominator in ( |S.15| ) . Instead of specifying the chiralities, we will consider 
positive as well as negative values of the momenta ps- Then we can write the fusion 
coefficients as a sum of 5-functions, e.g. 



NpiP2P3 ^Yl Hpi+P2+P3- (2n + l)po) (6.16) 

n=0 

Np^P2P3P4 + 1) S {pi + p2 + P3 + P4 - 2{m + l)po) , (6.17) 



m>0 



which we substitute in the general formula (|6.6| ). To write the final expression we must 
first evaluate the integral 



n-1 



= ^i^Pi - 2npo) X ^ {{pi +P2- {2k + l)pof + permutations) . (6.18) 



n>l i fc=0 

The expression multiplying the delta function simplifies to 

n-1 4 „ 4 



I](5Z(^'i+^'j)"^2/c + l)po) =n{Y,Pi~Pl) (6.19) 

k=0 3=2 i=l 
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where we used the identity (assuming J2Ps = 2npo) 



4 4 

$^(pi +Ps- npoY = J2p^s- nhl- (6.20) 

s=2 s=l 



Now we can write the exphcit expression for the r.h.s. of ( |6.6|) , which is again hnear in the 
momenta, 

j\^pi + ...+P4-2-2po ^ 

G{pi,P2,P^,Pa) = 7T~fY2 2^5(ptot - 2npo)x 

\ ) n=l 

1 . (6.21) 

n— 1 4 ^ ' 

f n . ^ —t \ 

n 



(2 + irl ^0) - 5Z 5Z l^'i + - + l)po\ 



k=0 s=2 



where ptot = Pi + P2 + Pa + P4 is the total momentum. Passing to world sheet CFT 
notations, we write ( |6.21|) at the critical points as 

CO / 4 \ n— 1 4 

G'(Pi, P2, P3, P4) ~ 5^ 5 5^ - 2neo [^Q - 5Z l^i + " + l)^ol 

n=l \ 1 / /c=0 s=2 

(6.22) 

This expression is identical with the expressions for the 4-point tachyon correlation function 
obtained for the "diagonal perturbation" of the non-rational Liouville gravity [|T|. More 
precisely ( |6.22| ) corresponds to the physical, symmetric in the four momenta correlators, 
extracted from the fixed chirality solutions of the functional relations. The intermediate 
states in ( |6.22D are spaced by 2eo and can be interpreted in terms of insertions of a screening 
vertex operator (tachyon) with matter component e^*'^"-^. This amounts to adding an 
additional "diagonal perturbation" to the Liouville term: 

^*=o = j d^a (A.e^^'^ + A^e^-o^e^^"^) . (6.23) 

In in the functional equations for the correlation functions were obtained by adding the 
action obtained by Liouville refiection, Aj^^^ = X^e^'^/^ + A^e^^'^^^e^'^. 

6.4- The four-point function of degenerate fields in the SRSOS model 

Here we assumed that pi = niipo are sufficiently small, so that the intermediate momenta 
are smaller than 1. The wave functions (|6.14 ), restricted to the discrete spectrum p G Pq'Zj+, 



^,„(,.) ^ ^ = ^'"(7°"^) . p = men, (6.24) 

sm(7rpoa;) 
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form closed Pasquier-Verlinde-like fusion algebra: 

^ {m-\-n— |m — 7T,| ) 

i^mpo'^npo = '^{m+n+l-2k)po ■ (6.25) 

k=l 

This is the algebra of the order operators. The multiplication rule in ( |6.25| ) coincides with 
the tensor product decomposition multiplicity of sl{2) irreps of dimension rrii 

|mi — m2| + 1 < rris < nii + m2 — 1 

and mi + m2 + ms = odd ' . (6.26) 

otherwise 

The 4-point multiplicity is correspondingly 

-^7711771277137714 ^ ^ -^771 1 777 2 -^77777737774 

m>0 (6.27) 
= I (min(mi + m2, ms + rn^) — max(|mi — m2|, |m3 — m4|) . 



The n-point multiplicities are given by the general formula ( p.20| ). The sum over x G Z+ 
can be replaced, up to an infinite volume factor, by an integral w.r. to the compact variable 
6 = irpox: 



N 

-''7771,..., 777^ 



/.27r 



To evaluate the sum in (|6.2| ), we first rewrite the expression as 

G(mi, m2, ms, m4) = ) ^ ^ (b.29) 



where the symbol )— /(m)— ( is defined similarly as in (|6.18|) 



• 771777,2 7713777 ^ ■ 



) f ij^^ ( • ^ ^ f ij^^ (-^77117772 777-^777777,37774 -^77717773777-^7712 7714 -^7771777,4777-^7] 

777 

(6.30) 



It is straightforward to check that 

4 



) 771 ( -^777 1 7772 77737714 ( ^ ^ -^777 1 777 2 777 3 7774 ) 



= 1 



4 



(6.31) 



= ^r77l777277737774 (Xl"^i -l). 



7 = 1 
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Inserting this in (|6.2| ) we obtain the foUowing simple expressions for the 4-point function 
of order operators. 



Gmim2m^m4 2 ~^ Po) -^7711771277137714 ^ ^ [-^77117712777 (^Po) -^77177737714 ~l~ perms.]^ 



771=1 

4 



2 

2 "^™'1"T'2 77137774 ( ^ ^~ -|^_£P0 PO ^ ^ ^7 ^~ P0-^771l777 2 77l377l4 j ' 

(6.32) 

In the two critical points, i = and £ ^ oo, the formula ( |6.32|) for the 4-point function 
can be compared with the expression obtained in [Q for the CFT with the "diagonal 
perturbation" ( |6.23| ) using the ground ring relations in the Coulomb gas approach. We 
find agreement, up to an overall constant. In the case of finite temperature, at the moment 
there are no world sheet CFT results to compare with. 

6.5. The 4-point function in ADE string theories 

At rational values of the background momentum, pq = 1/h, He the integers are 
restricted to the range 1 < nii < h — 1 and the s/(2) multiplicities ( |6.26| ), ( |6.28| ) are replaced 
by the Verlinde fusion multiplicities, i.e., the upper bound in ( |6.26|) becomes min(mi + 
m2 — 1, 2/i — mi — 7712 — !)• The corresponding relations ( |6.31| ) do not hold true anymore 
in general, however the initial quadratic formula (|6.2|), derived in |]13,B3|, survives and 



represents the tachyon 4-point function for pi = rriiPQ. This formula extends to the ADE 
cases. The 3-point multiplicities Nmim2m3 of the ADE series are real numbers (square 
roots of rationals), defining the structure constants of the Pasquier algebra [Q associated 
with any ADE graph, with the numbers corresponding to the set of exponents^. 

The partition functions of the ADE string theories on a torus are calculated in 
They are obtained by summing only over the order fields propagating along one of the 
circles of the torus. Nevertheless they reproduce the expressions obtained from the world 
sheet CFT [^, and there are no configurations that are not taken into account. This can 
be explained by the particular choice of the time slice, which goes along a domain wall. 
The torus is thus obtained by identifying the boundaries of a cylinder with Dirichlet type 
boundary conditions. 

The fact that the partition functions can be reconstructed only by taking into account 
the diagonal fields is also related to the fact that there is a world sheet CFT, the diagonal 
theory in |jl|] , in which these fields form closed algebra. One can give the following heuristic 
explanation of that. In the case of rigid geometry an order operator at the point r is 



10 



These constants coincide with the relative scalar OPE coefficients of the ADE theories [40|. 
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constructed by cutting out a circular domain having as a center this point and large 
compared to the lattice size. In general, there will be loops that enter and leave the 
domain, crossing its boundary. Therefore, when shrinking the domain to a point in the 
continuum limit, we will create also other operators as the thermal operators mentioned 
before. If we calculate a correlation function of two such operators, we will reproduce 
the standard OPE expansion of the diagonal fields (r = s) in CFT, which will include 
non-diagonal fields as well. Now consider the SRSOS model on a fluctuating lattice and 
repeat the procedure. In such a "liquid" lattice the symmetry is stronger and the notion of 
circular domain has no sense anymore. The domain surrounding the point r is formless; it 
has a single geometrical characteristics, its area and perimeter. We can restrict ourselves 
to the set of domains that do not cross any domain walls, i.e. with Dirichlet type boundary. 
Then the deflnition of the order operator by shrinking the boundary of the domain is such 
that in the 4-point function of order operators there will be only order operators in the 
intermediate channels. This is exactly what happens in the diagonal worldsheet string 
theory The deflnition of the order operators such that they form closed algebra is 
possible only in the theory coupled to gravity. 

7. Summary and discussion 

The theory we studied here - the SRSOS height model on random triangulations - is a 
good candidate for a generic microscopic realization of non-rational 2D gravity. At the 
rational values of the background momentum po it describes the A-series of the ADE string 
theories. The model is characterized by semi-inflnite target space, which we identifled with 
the Aoo Dynkin graph. The semi-inflniteness of the target space leads to identiflcation of 
the positive and negative momenta, which implies a reflection property of the correlation 
functions. It is also possible to construct a model that generalizes the height models of 
the D series coupled to gravity. The construction of the Doo height model is sketched in 
Appendix B. 

We formulated the SRSOS model as a conformal fleld theory on a Riemann surface 
representing an inflnite branched cover of the complex plane. The points of the Riemann 
surface label the FZZ branes in the SRSOS model. 

The calculation of the correlation functions is straightforward and it matches with 
the world sheet CFT predictions. Our results are compatible with the conjecture that the 
SRSOS model gives a microscopic realization of the "diagonal" CFT introduced in f^], in 
which the screening charges are the Liouville dressed identity operator (with momentum 
Po) and its charge conjugated (with momentum —po). 

The SRSOS model gives formulas only in the physical regions of momenta, for which 
the tachyons describe local observables on the world sheet. The expression is given by 
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a sum of terms associated with intermediate channels. Each such term is a polynomial 
of the momenta of the particles involved. These polynomials are higher order than the 



polynomials that occur in the iteration procedure of Di Francesco and Kutasov [|T5[. In 
the case of non-rational theory there are relations that allow to lower the degree of the 
polynomials. Using them we have transformed the discrete model 4-point functions into 
expressions comparable with the CFT solutions of the ring relations. This property does 
not seem to extend to the case of the minimal matter theory with rational 6^ - a case which 
is especially difficult to be analyzed in the CFT approach. Thus we expect that the 4- 
point functions in this case are given by expressions similar to the initial quadratic formula 
( |6.2|) . The discrete approach provides as well the analogs of the ADE correlators, which 
cannot be treated in the CFT approach, as the bulk ground ring relations, expressing the 
fundamental isospin 1/2 fusion relations, do not apply directly to the non-diagonal D,E 
cases. 

In this paper we restricted ourselves to the correlation functions of order operators. 
In fact our general formula ( |6.32|) is the generating function of the correlators of 4 order 
operators and an arbitrary number of Liouville dressed thermal operators $1,3, which are 
given by the expansion coefficients of the thermal coupling t. The meaning of the results 
for finite temperature will be discussed in detail elsewhere [E5 . 
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Appendix A. Loop gas representation of the correlation functions 
A.l. Correlation functions of order operators in the SOS model 

The order, or vertex, operator in the SOS model is defined as the effect of inserting the 
wave function 

^p{x) = e^^(Po-P)^ = e^*'^^ (A.l) 

at some point of the planar graph. We will refer to p as a target space momentum and to 
9=1 {Po ~ P) 3-3 electric charge. The correlation functions of the vertex operators (|A.1| ) 
can be formulated as the partition function of the loop gas with some of the loop weights 



40 



modified. The rule is that the weight of a loop encircling a total charge q^'^ = ^{po — 
is equal to 

2 cos( = 2 cos[7r(po - 2?'°)] . (A.2) 

Since the reference point (the infinite point for the plane) can be put anywhere, this rule 
is consistent only if it gives the same result with q^^ replaced with the total charge 
outside the loop. This is guaranteed by the charge conservation condition q^^ + q°'^^ — po 
or, in terms of momenta, p^^ + p°^^ = 0. 

• One-point function 

There is only one vertex operator, •i/'_pg(x) = e^*'^^"^, with non-trivial one-point function. 
This operator is related to the identity operator by the charge refiection 

Q ^Po-Q or p^ -p- (A.3) 

The charge of this operator compensates the total background charge —po for the sphere. 
The correlation function {'4>-po) is equal to the partition function of loop gas with weight 
2 cos(7rpo) per loop. 

• Two-point function 

The two points can be connected by a line that intersects each loop at most once. The 
two-point function ^ ■^p'^_p ^ is obtained by changing the weight of the non-contractible 
loops, i.e., those that intersect the line, to 2cos7rp. 

• n-point function {n > 3) 
Consider the n-point function 

Gp^,...,p^ = ^ •0pi'0P2 • • •V'pn ) , ^Pi = {n~ 2)po modulo 2 (A.4) 

i 

where pi, . . . ,Pn-i are positive and Pn is negative ("chirality rule"). The n points where 
the operators are inserted can be connected by a set of oriented lines that form a tree, with 
the condition that the tree intersects each loop at most once. An example for n = 4 is 
given in Fig. 10, where the fourth point is at infinity. To each line of the tree we associate 
a charge. The charges associated with the external lines are those introduced by the order 
operators. The charges associated with the internal lines are determined by the charge 
neutrality of the vertices of the tree. In terms of momenta the neutrality condition for a 
vertex with k lines is written as J2iiPi ~ Po) + 2po = modulo 2. Then the loops crossed 
by a line I of the tree change from 2cos(7rpo) to 2cos(7rp/), where pi = po — 2qi is the 
momentum associated with this line. The above argument hold unchanged if instead of 
point-like insertions one considers finite boundaries. 
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Fig. 10: A loop configuration in presence of four operator insertions. 



A. 2. Correlation functions of order operators in the SRSOS model 

Unlike the SOS model with po ^ 0, the SRSOS model has a well defined partition function 
on the sphere without boundaries. Summing over all height configurations yields a factor 
2 cos Tcpo per loop. Thus the partition function of the SRSOS model on the sphere coincides 
with the one for the 0{n) model with n = 2cos(7rpo)- 

The order operators are introduced as insertions of the wave functions 

, . , sm(Tcpx) ^ , . ^, 

sm^irpox) 

If we restrict the spectrum of momenta to 

p = mpo, m = l,2, ... (A. 6) 

then the functions ([A.5| ) generate the su{2) character ring. 

• One-point function 

The only non-trivial one-point function is that of the identity operator, p = po- One can 
perform the sum of the heights by repeatedly applying the relation starting with 

the domains that do not contain loops. As a result, the expectation value (i^po) is equal 
to the partition function of the loop gas on a sphere with one marked point, with weight 
2cos(7rpo) for each loop. 

• Two-point function 

The insertion of two operators (|A.5| ) at two points changes the weights of the non- 



contractible loops to 2cos7rp. The weights of these loops follow from the identity ( 2.11J ) 
with Po replaced by p. This two-point function is defined for any real value of p. 

• n-point functions (n>3) 
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The n-point functions in the SRSOS model are weU defined when the momenta belong to 
the spectrum ( |A.6|) , 



C^i.-.-.m™ = (V'piV'pa ••■V'p„), Pi = miPQ. (A.7) 

One can again describe the weights by drawing a tree connecting the n points and such that 
each loop is intersected at most once. The momenta mi, ...^rriji in (|A.7D are associated 



with the external lines of the tree. Unlike the SOS case, here the momenta associated 
with the internal branches are not determined uniquely. When calculating the correlation 
function with n > 3, a sum over these momenta is to be performed. The weights in 
the momentum space are associated with the vertices of the tree. If m^^, ...,mi^ are the 
momenta associated with a /c-vertex, then the weight of this vertex is given by the su{2) 
multiplicity Nmi^,...,mi^j which is the Fourier image of the weight ( |2.9| ) with n = k. 

It follows from the explicit form of the su{2) multiplicities that if the n external 
momenta satisfy the neutrality condition 

mi + . . . + mn-i - mn = {n - 2), (A.8) 

then the internal momenta are uniquely determined and the n-point correlation function 
( |A.7| ) coincides with the n-point function in the SOS model with n — 1 positive momenta, 



Pi = miPo (1 < k < n — 1) and one negative momentum, p^ = —m^Po- Thus the n- 
point functions in the SRSOS model satisfying ( |A.8| ) are identical to the corresponding 



amplitudes in the SOS model that satisfy the "chirality rule" . 
In the general case eqn. ( |A.8|) is replaced by 



mi + . . .-\- mn-i - mn = {n-2) + 2k, keZ+. (A.9) 

The contribution of given loop configuration to the n-point function can be written as a 
sum of SOS n-point functions with k extra operators i/j-p^ix) — 6*^^°^, associated with 
the domains containing the vertices of the tree. 



Appendix B. Expression for the operators a[ ^ in terms of the collective field 

We will use the expression ( |5T^ ) of the creation operators aj^^ as contour integrals of 
the effective potential T^iz). The operators a]^^ is proportional to the coefficients in 
the expansion of the effective potential in the half-integer powers of z + M. In the r 
parametrization 

yiT^=v^coshf . 

Instead of expanding Tj:{z) at ^ = —M we will evaluate the integrals ( p.2|) by deforming 
the contour C_ to the contour C-(- going around the second cut of T^lz) along the interval 
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[M, oo] as in Fig. 5b. Deforming the contour C_ to C+ (with clockwise orientation) we 
get, up to an infinite constant, 



1 



T—i 



s. 



4,1^) = f ir- 



9: 



dz 



dz 



fB.ll 



The denominator in the last integral is regular and we express the discontinuity of the 
field V^{t{z)) in terms of $(r(z)). To evaluate the integral we go to the r parametrization 



1 



oo „ 
\ - U 



V29.. 



s n=0 ^' 



dr dr sin rrxdr cos Tvdr $ (r) 
TT Vl + coshr — 2u 



(B.2) 



In this step we have retained the negative frequency part of the field, i.e., ( [B.2D is un- 
derstood to hold applied on the right vacuum |0oo)- The projections ([B.2|) are seemingly 
non-local, but they are actually given by polynomials of the derivatives of the holomorphic 
field (/?(t) at r = ±inx. To see that we evaluate the integral for each Fourier mode in the 
expansion ( [4.16|) of ip~ (r) using the integral formula for the Legendre function of first kind 

m- 

dr coshi^r P_i+,(l-2w) 

(B.dj 



TT Vcosh T + 1 — 2u 

■1 



V2 



cos TTZ^ 



k=0 



k\ 



k\ 



1,2,... 



The integral over v amounts in replacing u dr and setting r = 0, i.e., we get 



2 ^ TT^ OrSm TrxOr ^' (tj 



T = 



;B.4) 



(B.5) 



or, transforming to the momentum space and using the shorthand notation (|4.8| ) for the 
differential operator, 

{aiip)-5k,ig:'5^'\p-po))\Ooo)=9sFk{dr)drip-ip.r) J0,J 



T = 



sinhr^r {r + in.p) \0oo) ■ 



(B.6) 



T = 



On the other hand assuming that we can identify |0oo) with the right twisted vacuum 
l^tw)' can replace in the r.h.s. dr~ (r + iiv) with the twisted negative mode part. Using 
that the decomposition ( |5.2| ) takes the form 

1 . . S... (-2)"-i(sinh. 



art(r + 



in 



n>0 



r Sx ^ 
9s 



T\2n-1 
2' 



(2n- 1)!! 



(B.7) 



and plugging (|B.7|) in the identity ( [B.6|) we see that it implies that the operators Fk act 
as projectors for r = 0, namely. 



Fk{dr) sinhr 



2"-i(_l)'^(sinh|)2"-i 



|r=0 



Sk, 



n > 



fB. 



(2n- 1)!! 

This property is also checked independently to hold true, which in turn ensures the iden- 
tification of the right vacua above. 
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Appendix C. The Doo model 

The target space of the D^o model is labeled by the positive integers x G Z+ and the 
two extremities of the 'fork' a; = and x — ^ (Fig- H)- The adjacency matrix is given by 

= (^x.x'+i + x,x' en, 

(CI) 




1 2 3 4 ... 



Fig. 11: The target space of the Doo model 

(p) 

Its eigenvectors Sx \ Q <p <'k, are given by 

^(,) ^ r cos(7rpx) ifxGZ+; q < p < ^. (C.2) 

^ [ ^ itx = Oorx = 0. — — \ J 

The zeroth eigenvalue = or p = | is twice degenerated. The second eigenvector, 
which we denote by 5"^^/^^ has components 

[ \ ifx = 0; 
S^^'^^=l -i ifx = 0; (C.3) 
[ if X G Z+. 

For po = 2n+2 target space can be restricted to the Dn Dynkin graph. 
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